Preprint typeset in JHEP style - HYPER VERSION 



Modeling Heavy Ion Collisions in AdS/CFT 



00 ■ 

o : 

Javier L. Albacete, Yuri V. Kovchegov, Anastasios Taliotis 

Department of Physics, The Ohio State University, Columbus, OH 43210, USA 



E-mail addresses: albacete@mps.ohio-state.edu, yuri@mps.ohio-state.edu 



taliotis . l@osu.edu 



Abstract: We construct a model of high energy heavy ion collisions as two ultrarelativistic 
shock waves colliding in AdSs. We point out that shock waves corresponding to physical energy- 
momentum tensors of the nuclei completely stop almost immediately after the collision in AdSs, 
which, on the field theory side, corresponds to complete nuclear stopping due to strong coupling 
effects, likely leading to Landau hydrodynamics. Since in real-life heavy ion collisions the large 
Bjorken x part of nuclear wave functions continues to move along the light cone trajectories of the 
incoming nuclei leaving the small- a; partons behind, we conclude that a pure large coupling approach 
is not likely to adequately model nuclear collisions. We show that to account for small-coupling 
effects one can model the colliding nuclei by two (unphysical) ultrarelativistic shock waves with zero 
net energy each (but with non-zero energy density). We use this model to study the energy density 
of the strongly-coupled matter created immediately after the collision. We argue that expansion 
of the energy density e in the powers of proper time r squared corresponds on the gravity side to 
a perturbative expansion of the metric in graviton exchanges. Using such expansion we reproduce 
our earlier result [1] that the energy density of produced matter at mid-rapidity starts out as a 
constant (of time) in heavy ion collisions at large coupling. 
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1. Introduction 

There is a mounting phenomenological evidence coming from RHIC data for a strongly coupled 
medium created in heavy ion collisions [2-19]. Ideal hydrodynamics simulations have been extremely 
successful in describing data generated in heavy ion collisions at RHIC [2-9] . These analyses require 
very small shear viscosity [10], indicating that the medium is strongly coupled [11-15], and a very 
short thermalization time of the initially produced system, of the order of 0.3 4- 0.6 fm/c [2-9]. At 
the same time, many bulk features of heavy ion collisions at RHIC which are sensitive to the initial- 
time dynamics, such as the energy, rapidity and centrality dependence of particle production are 
very well-described in the weakly-coupled framework of the Color Glass Condensate (CGC) [20-38] 
(for a review of CGC see [39-41]). In the CGC approach a heavy ion collision releases the small 
Bjorken-x partons in the nuclear wave functions, which quickly go on mass shell and become real: 
the resulting particle distribution in highly anisotropic in momentum space [42,43]. A question 
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then arises about how this initially weakly-coupled highly anisotropic system becomes isotropic and 
equilibrates very quickly becoming a strongly-coupled, possibly thermal, medium. 

So far, conventional perturbative descriptions of thermalization of the produced medium [44-47] 
have not been able to account for the short thermalization time of the order of a fraction of a fermi 
required by hydrodynamic simulations to describe the data. Perturbative thermalization in heavy 
ion collisions was also studied in [43,48,49], where it was concluded that perturbative thermalization 
scenarios are not likely to be compatible with RHIC data. This leads one to conclude that it 
is highly likely that non-perturbative large-coupling QCD effects are responsible for the apparent 
thermalization observed in RHIC data. While the research on perturbative thermalization scenarios 
along the lines outlined in [44-47] is vigorously pursued in the community, we believe it probable 
that the dynamics of the medium produced in a heavy ion collision at RHIC proceeds as follows. 
The medium starts out being weakly-coupled, being well-described by CGC. After a very short 
proper time, of the order of r> 1/Q S 0.2 4- 0.3 fm/c (and possibly at r ps* 1/A qcd ps 1 fm/c) 
the coupling becomes strong. Strong coupling effects are likely to quickly thermalize the medium, 
allowing for its hydrodynamic description. 

Unfortunately to date there is no consistent way to describe both the weakly-coupled and the 
strongly coupled dynamics of QCD medium in a unified framework. However, as we argued above, 
to understand the general physical nature of thermalization (and, more importantly, isotropization) 
of the medium, and to see whether isotropization and thermalization take place at all, a purely 
strong coupling approach seems, a priori, appropriate. Indeed at strong coupling QCD becomes non- 
perturbative and no controllable dynamical calculation appears to be possible. Instead one could use 
the Anti-de Sitter space/conformal field theory (AdS/CFT) correspondence [50-53] to understand 
the same process for Af = 4 super- Yang-Mills theory. This correspondence, and in particular the 
gauge-gravity duality which follows from it, allows one to understand strong coupling effects in 
such QCD-like theories as Af = 4 super- Yang-Mills theory using super-gravity in 5 dimensions. 
AdS/CFT correspondence has been useful in providing insight in the behavior of the shear viscosity 
in strongly-coupled gauge theories [11-15] along with providing other interesting results on the 
evolution of the medium created in heavy ion collisions [1,54-67]. 

The goal of this paper is to make progress in constructing a dual geometry in AdSs space for 
heavy ion collisions with the goal of understanding the onset of isotropization and thermalization 
of the produced medium. The previous paper on the subject by two of the authors [1] studied the 
very early time dynamics of the medium produced in the collisions. It was shown that, assuming 
rapidity-independence of the produced medium and assuming non-negativity of its energy density, 
one would obtain that the energy density of the produced medium should start out as a constant of 
time at very early times immediately after the collisions. This implied that the longitudinal pressure 
of this early-time medium is negative and the medium is thus highly anisotropic. This behavior 
is similar to that of the weakly-coupled CGC medium at early times [32,68,69]. The problem of 
isotropization and the onset of Bjorken hydrodynamics [70] in this framework can be formulated 
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as the question about understanding the transition from the negative longitudinal pressure of the 
medium at early times to the positive longitudinal pressure (comparable to the transverse pressure) 
at late times [43]. In [54,56,71] it was shown that the dynamics of a strongly-coupled rapidity- 
independent medium leads to Bjorken hydrodynamics behavior at late proper times: it is therefore 
likely that isotropization transition takes place at some intermediate time between the early-time 
dynamics of [1] and the late time dynamics of [54,56,71]. 

To better understand this transition one needs to find the energy-momentum tensor of the 
medium at a later times than considered in [1]. Unfortunately, a consistent expansion of the 
energy-momentum tensor in the powers of proper time r requires some knowledge of the geometry 
dual to the colliding nuclei. (In [1] nothing was assumed about the colliding nuclei, except that 
they lead to rapidity- independent distribution of matter.) Thus in this paper we try to construct 
a geometry dual to the collision of the two nuclei. We model two nuclei as shock waves in AdSs. 
Modeling nuclei with shock waves has previously been considered in [66] for AdS3 corresponding to 
gauge theory in two space-time dimensions and in [61,67] for AdSs. 

The paper is organized as follows. We begin by spelling out some general formulas used in the 
paper in Sect. ||. We proceed by setting up the problem in Sect. |[ We start with two colliding 
shock waves, each given by the metric like that shown in Eq. (|3.1|) (see [54]). We argue that this 
metric corresponds to a single graviton field, as shown in Fig. ^|. We then argue that a consistent 
expansion of the metric at the time after the collision can be constructed by considering higher 
order graviton exchanges between the boundary and the bulk, as shown in Fig. We construct a 
general next-to-leading order perturbative contribution to the metric in graviton exchanges in Sect. 
| The solution is given by Eqs. Q, (TO , ( CT , CT , ( TO , ( H ), (EH), ( CT , ( PI - 

Using the obtained solution we study the collision of two physical shock waves in Sect. |5]. We 
conclude that the shock waves, and the nuclei in the boundary theory, completely stop very shortly 
after the collision, after a time of the order of inverse typical transverse momentum scale in the 
problem, as given in Eq. ( [5.19| ). We interpret this result as creation of a black hole, similar to what is 
suggested for collisions of particles at transplanckian energies [72-74], though our black hole would 
be created in the bulk. On the gauge theory side this implies that strong coupling effects would 
completely stop the nuclei shortly after the collision, on the time scale less than or equal to 1 fm/c. 
Such strong coupling effects are likely to thermalize the system soon after the stopping, leading to 
Landau hydrodynamic description of the system [75]. This thermalization scenario is very different 
from the onset of Bjorken hydrodynamics outlined above. It is possible in principle and has been 
advocated in the literature [76]. 

However, we believe that the Landau hydrodynamic scenario is not likely to be relevant for heavy 
ion collisions. This claim is supported by the following observations. On the one hand, the ideal 
Bjorken hydrodynamics has been extremely successful in describing the particle spectra and elliptic 
flow [2-9] . On the other hand, the agreement with experimental data of the predictions based on the 
weakly-coupled CGC approaches to description of particle multiplicities in nuclear collisions [33-35] 
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indicates that the complete nuclear stopping does not happen in the actual collisions. Indeed, 
in high-energy scattering at small coupling the hard (large Bjorken-x) parts of the nuclear wave 
functions simply go through each other without recoil. In the quasi-classical CGC limit this leads 
to rapidity-independent, Bjorken-like picture of particle production in heavy ion collisions [26,27, 
77-79]. Thus it appears that weak coupling effects are a key ingredient for a proper description 
of the space-time structure of heavy ion collisions, even in the medium becomes strongly coupled 
shortly after the collision. 

As we do not know how to model weak coupling effects in AdS/CFT, we try in Sect. || to mimic 
them by introducing unphysical shock waves with non-positive-definite energy density. Such shock 
waves are indeed unphysical and can not follow from the underlying string theory. Even for the 
gravity in the bulk these shock waves can only serve as sources external to the theory. An example 
of the energy- momentum tensor of such shock wave is given in Eq. ( |6.3| ). Using the general solution 
of Sect. [|, we construct the early time energy density and pressures of the medium produced in the 
"collision" of two of such shock waves, which are shown in Eq. ( |6.7| ). As can be seen from Eq. ( |6.7| ), 
the energy density starts out as a constant in time. We have thus reproduced the results of [1], 
but now in a more dynamical setting. The problem of thermalization formulated with the help of 
these unphysical shock waves is probably that of isotropization of the produced medium to achieve 
Bjorken hydrodynamics [70], as described above and considered in [1,43]. 

We discuss possible higher order corrections to the result of Eq. ( |6.7| ) in Sect. and observe 
that a dilaton field may need to be introduced at higher orders to account for initial non-equilibrium 
between chromo-electric and chromo- magnetic modes in the medium (see Eq. ( [7.6[ )). We conclude 
in Sect. § by restating our main results. 

Our solution found in Sect. ^ is general and is valid for any shock wave profile, unlike the 
solution found in [67] which works only for delta-function shock waves. In the particular case of the 
delta- function shock waves, our solution reduces to that found in [67]. The general nature of our 
solution allowed us in Sect. |5] to reach a more general physical conclusion that the collision of any 
two physical shock waves (with positive definite energy density) leads to complete stopping of the 
shock waves after the collision, probably leading to the formation of a black hole. In particular this 
conclusion applies to the delta-function shock waves used in [67]. Our solution also allows us to 
tackle unphysical shock waves in Sect. |5], which are impossible to handle in the formalism of [67]. 
We should also note that the stopping of physical shock waves found here does not happen in AdS 3 
(see [66]), since in 1+1 dimensional gauge theory the nuclei are point particles and stopping for 
them is impossible, as there is no transverse directions in which the momentum could be channeled. 

2. Some Generalities 

Throughout this paper we will work with the metric of AdS5 written in terms of Fefferman-Graham 
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coordinates [80] 



7-2 

ds 2 = — {g^(x,z)dx»dx u + dz 2 } (2.1) 



where /i, v run from to 3 and z is the coordinate describing the 5th dimension. The boundary of 
the AdS space is at z = and L is the curvature radius of the AdS space. 

According to holographic renormalization [81], if one expands the 4-dimensional metric g^ u {x, z) 
near the boundary of the AdS space 

~9^z) = ~g^(x)+z 2 ~g^(x)+z i ~g^(x) + ... 1 (2.2) 

then, for Minkowski metric g$ (x) = r]^, one gets gffl(x) = and the expectation value of the 
energy-momentum tensor of the gauge theory is 

N 2 

Below we will use the light cone coordinates 

/Y»0 — I — ,-y»3 

_i_ X —I— Jj . . 

x = — 

where x 3 is the collision axis of the two colliding nuclei. In these coordinates the empty AdSs metric 
is 

L 2 

ds 2 = — { -2 dx + dx~ + dx\ + dz 2 } , (2.5) 

z 

where dx\ = (dx 1 ) 2 + (dx 2 ) 2 with x 1 and x 2 the transverse dimensions which we will denote using 
Latin indices, e.g. x l . To describe nuclear collisions we will also use the proper time 



and space-time rapidity 



Einstein equations in AdSs are 



V2 x+x- (2.6) 



R^-^g^R + A c g^ = (2.8) 

where g^ is the full 5-dimensional metric of the AdSs space, R is the scalar curvature and A c is 
the cosmo logical constant. For AdSs we have 

Ac = ~ (2.9) 
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and Eq. (j2.8|) gives 



R = - 



20 
1?' 



(2.10) 



Eqs. (pD, yield 



Ry,u "I - ^2 SVw 0" 



(2.11) 



Later in the paper we will also discuss the dynamics of a dilaton field if coupled to gravity in 
AdS 5 . In the presence of dilaton Eq. ( |2.8| ) is modified to (see e.g. [71,82,83]) 



1 1 1 

Riw ~ ^9^uR + A c = - df.ip d u <f - - g^ d p ip <9V 



Eq. ( |2.12|) can be simplified to give 

4 1 

Rnu + J2 g^u = 2 d n<P dvV- 

The dilaton itself obeys the Klein-Gordon equation in curved space-time 



(2.12) 



(2.13) 



(2.14) 



3. Setting up the Problem 

Out goal is to construct a metric in AdSs which is dual to an ultrarelativistic heavy ion collision 
as pictured in Fig. |I[ Throughout the discussion we will use Bjorken approximation of the nuclei 
having an infinite transverse extent [70] and being homogeneous (on the average) in the transverse 
direction, such that nothing in our problem would depend on the transverse coordinates x 1 , x 2 . 

We start with a metric for a single ultrarelativistic nucleus moving along a light cone. As was 
noted by Janik and Peschanski [54] the following metric gives a solution of Einstein equations in 
AdSs in Fefferman- Graham coordinates [80] 

ds 2 = — j-2dx + dx~ + — (T__(aT)) z 4 dx~ 2 + dx\ + dz 2 j . (3.1) 



Eq. ( p.l| ) is a solution of Einstein equations ( |2.11| ) for any expectation value of the energy-momentum 



tensor of the nucleus in four dimensions (T {x~)), as long as it is a function of x~ only [54]. The 

factor of 2tt 2 /N 2 is due to Newton's constant [81]. For an ultrarelativistic nucleus with infinite 
transverse extent moving along the x + axis (see Fig. P the leading components of the energy 
momentum tensor depend only on x~ . Hence the metric in Eq. ( |3.1|) adequately describes such a 
nucleus, though does not restrict the dependence of (T (x~)) on x~ . 
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nucleus 1 



nucleus 2 



Figure 1: The space-time picture of the ultrarelativistic heavy ion collision in the center-of-mass frame. 
The collision axis is labeled x 3 , the time is x . 

While the metric (|3.1| ) is an exact solution of the non-linear Einstein equations ( |2.11| ), it can also 
be represented perturbatively as a single graviton exchange between the source nucleus at the AdS 
boundary and the location in the bulk where we measure the metric/graviton field. This is shown 
in Fig. 0, where the solid line represents the nucleus and the wavy line is the graviton propagator. 
Incidentally a single graviton exchange, while being a first-order perturbation of the empty AdS 
space, is also an exact solution of Einstein equations. This means higher order tree-level graviton 
diagrams are zero. It is interesting to note that a similar property has been observed for gauge 
theories in covariant gauge [24,84], where the exact solution of Yang-Mills equations with a single 
ultrarelativistic nucleus as a source is given by a single gluon exchange. 

As one can see comparing the metric ( |3.1|) with the diagram in Fig. [|, each graviton-nucleus 
vertex gives a factor 

~ ±(T-{x-)). (3.2) 

c 

If the nuclear energy-momentum tensor is iV c -independent, then in the large- N c limit the factor in 
Eq. (|3.2p would be small and one could envision perturbative expansion in this parameter for the 
problem of collision of two nuclei. On the other hand, gauge-gravity duality is valid only in the 
large- N c limit: hence we need (T (x~)) ~ N% to avoid having iV c -suppression for higher-order 
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N 



2 <T__> 



graviton exchanges. This could be achieved by imagining a nucleus with nucleons made out of 

N 2 — 1 "valence" gluons each. Then (T (x~)) ~ N 2 and multiple graviton exchanges will not be 

iV c -suppressed. 

However, one may then worry that the expansion parameter 
also ceases to be small. Nevertheless it makes sense to expand in 

powers of (T (%~)), as usually (T contains some momen- 

turn scale characterizing the nucleus such that one can keep track 
of the powers of this scale. Hence even if the expansion parameter 
is not small, the expansion is still well-defined and can be kept 
track of. 

For instance, in the original proposal of Janik and Peschanski 
[54], the energy-momentum tensor due to valence quarks in the 



ultrarelativistic nucleus was taken to be 

(T__0O) = ^N 2 J{x-) 



(3.3) 



with n a scale having dimensions of mass cubed. Starting with 
the energy-momentum tensor of a single ultrarelativistic particle 
and performing the averaging along the lines shown in Appendix 
[A| gives 

/I oc p + A 2 A 1/3 



Figure 2: A representation of 
the metric ( |3.1|) as a graviton 
(wavy line) exchange between the 
nucleus at the boundary of AdS 
space (the solid line) and the 
point in the bulk where the met- 
ric is measured (denoted by a 
cross) . 

(3.4) 



where p + is the large light-cone momentum of the nucleons in the nucleus, A is the atomic number 

and A is some transverse momentum scale. In this case the expansion in powers of (T (x~)) 

translates into the expansion in the powers of /j,, which can be systematically resummed (see e.g. 
[67]). 

Alternatively one could argue that at strong coupling the energy-momentum tensor is dominated 
not by valence quarks, but by the strong gluon fields of the nucleus. One can argue, based on 
conformal invariance, that the coordinate dependence of the energy-momentum tensor of such a 
strong gluon field in Af = 4 SYM theory is the same as that of weakly coupled electromagnetic 
fields [85-87]. Performing a classical electrodynamics calculation of the energy-momentum tensor 
of a point charge, averaging over all transverse coordinates and summing over all nucleons yields 
(see Appendix [A| for details) 



(T__(aT)) = V\K 2 N 2 c A l /H(x- 2 ) 

with A some transverse momentum scale. At weak coupling in classical electrodynamics 
(T (x~~)) ~ g 2 ~ A with A the 't Hooft coupling 



(3.5) 



A 



g 2 N c . 



(3.6) 
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Guided by the calculation of the heavy quark-antiquark potential at strong coupling in [88] we have 
replaced the coupling constant A by v^A in Eq. ( |3.5| ). However, the exact power of A in Eq. (|3.5| ) 
does not alter the subsequent discussion. 

In case of the energy-momentum tensor in Eq. (|3.5|) one can construct an expansion in the 
powers of A 2 , which is again well-defined. 




• • • 



~ N 2 <t; ++ > ~ N 2 <t; ++ > 

ABC 



Figure 3: Diagrammatic representation of the metric in Eq. (3.7). Wavy lines are graviton propagators 
between the boundary of the AdS space and the bulk. Graphs A and B correspond to the metrics of 
the first and the second nucleus correspondingly. Diagram C is an example of the higher order graviton 
exchange corrections. We calculate the contribution of this diagram below in Sect. ||. 



Using the perturbative expansion in the powers of the energy- momentum tensor, one can con- 
struct the metric dual to a heavy ion collision. At the lowest non-trivial order we begin by writing 
the metric as 

L 2 ( 2-K 2 2 7T 2 

ds 2 = —< -2dx + dx~ + dx 2 ± + dz 2 +—— (T^-tx-)) z 4 dx~ 2 + —— (T 2++ (x + )) z 4 dx +2 
z 2 i ~ N 2 ' N 2 

+ higher order graviton exchanges j> (3-7) 

where (Tj and (T 2++ (x + )) are the energy-momentum tensors of the colliding nuclei 1 and 2 

as shown in Fig. [I]. The metric in Eq. ( |3.7| ) is that of two colliding shock waves in AdSs. Higher order 
graviton exchanges will modify the shock waves after the collision and generate energy-momentum 
tensor of the matter produced by the collision in the forward light cone. The metric of Eq. (|3.7| ) 
is our formulation of the problem of heavy ion collisions in AdS. Similar metrics were previously 
considered in modeling heavy ion collisions in AdS3 in [66] and in AdSs in [61,67]. 
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The metric of Eq. (|3.7|) is illustrated in Fig. |^. The first two terms in Fig. ^| (diagrams A and B) 
correspond to one-graviton exchanges which constitute the individual metrics of each of the nuclei, 
as shown in Eq. ( |3.1| ). Our goal below is to calculate the next order correction to these terms, 
which is shown in the diagram C in Fig. ^. Higher order graviton exchanges would necessarily 
involve both nuclei: as the metric ( |3.1| ) is an exact solution of Einstein equations, all higher order 
graviton exchanges with a single nucleus are zero. Indeed, solving Einstein equations order-by-order 
in the graviton exchanges one could reconstruct any higher order term in the series of Fig. |||. In 
the calculations below we will restrict ourselves to diagram C in Fig. |3|, which is the first correction 
to the sum of the metrics of the two nuclei, and leave calculation of the higher orders for future 
projects. 



4. General Perturbative Solution 



Here we will calculate the diagram C in Fig. |^ by solving Einstein equations ( 2.11] ) perturbatively. 



Define normalized light-cone components of the nuclear energy-momentum tensors by 

ti(x-) = %{Ti-{x-)) (4.1) 

* c 



and 



t2{x + ) = ^(T 2++ (z+)>. (4.2) 

c 



Using these definitions we rewrite the metric in Eq. ( |3.7|) as 

ds 2 = ^- j -2dx + dx~ + dx\ + dz 2 + t l (x~)z i dx~ 2 + t 2 {x + )z A dx +2 + o(tit 2 )^. (4.3) 

Notice that there is no higher order corrections containing only powers of t\(x~) or of t2{x + ): they 
are zero since the single nucleus metric ( |3.1|) and its analogue for nucleus 2 are exact solutions of 
Einstein equations. 

We denote by g^ v the metric in AdS 5 space dual to heavy ion collisions that we are trying to 
construct. Then Eq. (|4.3j ) can be written as 



ds 2 = g iiV dx^dx v (4.4) 

with /i, v running from to 4. The order-by-order perturbative solution of Einstein equations is 
obtained by expanding the metric around the empty AdS 5 space 

9,u = ff{2 + ff{5 + $ + ■■■ ( 4 -5) 
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where the metric gffl corresponds to n graviton exchanges. For the energy-momentum tensor of 
Eq. (|3.4j) the series in Eq. ( |4.5|) corresponds to expansion in the powers of fi, while for the energy- 
momentum tensor of Eq. ( |3.5| ) the series in Eq. ( }4.5| ) is an expansion in powers of A 2 . 

Here gfij is the metric of the empty AdSs space with non-zero components 



(°) = (°) = -— (°) = §..— i j = 1 2 (°) = — (4 6) 



g$ is the first perturbation of the empty AdSs space due to the two nuclei 

g ( l = h{x-) z\ g+i = t 2 (x + ) L 2 z 2 (4.7) 
with all the other components zero. 

(2) 

We want to find the next non-trivial correction g^J. By the choice of Fefferman- Graham 
coordinates one has g Zfl = g^ z = exactly for \i ^ z and g zz = L 2 /z 2 . Hence the non-trivial 

(2) 

components of g)J are those for //, v = 0, . . . , 3. Due to translational and rotational invariance of 

(2) 

the nuclei in the transverse direction g\- ~ dij. We thus parametrize the unknown components of 

(2) 

g^v as 



J2 j2 

9-- = ~t f(x + ,x~,z), gfl = — f(x + ,x~,z), 
z z z z 



1 L 2 L 2 

9+- = ~2 9(x + , x~,z), g-f = h(x + ,x~, z) 5 i} (4.8) 



with /, /, g and h some unknown functions. Imposing causality we require that functions /, /, 
g and h are zero before the collision, i.e., that before the collision the metric is given only by the 
empty AdS space and by the contributions of the two nuclei Eq. (|4.7|) . Also, according to general 
properties of g^ v outlined in Sect. |2] (see [81]), we demand that /, /, g and h go to zero as when 
z -> 0. 



Using Eqs. (|4.6|) , (|4.7f) , and ( |4.8|) in Eq. (fig ), plugging the latter into Einstein equations 
( [2. 11| ) and keeping only the terms up to and including the order gff we obtain the following set of 



equations for /, /, g and h labeled by the Einstein equations components: 
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-+ 



[zz) 

:+z) 



3 _ ^ _g 

t/ Z J Z* Z XX } 

fz~ -fzz - h x + x + = 0, 



2^ " 2 
5 1 



1 1 

, , i . h-z , 9zz „ X~ <^X+ 2 

4 z z 4 2 
4z 6 t 1 (x-)t 2 (a; + )-iz 8 t / 1 (x-)4(a; + ), 



# 2 + 5/i z - zh zz + 2zh x+x - = &z 7 ti{x )t 2 (x + ), 
g z + 2h z - zg zz -2zh zz = -32 z 1 t x (x~) t 2 (x + ). 
_1 f 1 _ _ 

2 «/ X~^~ Z ii Qx Z X 

- ~ fx- z - ^9x+ z - K+ z = -z 1 ti(x-)t' 2 (x + ). 



-z 1 t' 1 (x-)t 2 (x- 



- l -~f 



(4.9a) 
(4.9b) 

(4.9c) 

(4.9d) 
(4.9e) 

(4.9f) 
(4.9g) 



The subscripts z, x + and x~ indicate partial derivatives with respect to these variables. 

To solve Eqs. (|4.9|) begin by solving Eq. (|4.9d|) for g z and substituting the result into Eq. (|4.9e|) . 
This gives 

— 3 h z + 3 z h zz — z 2 h zzz + 2 z 2 h x + x - z = 16 z 1 ti(x~) t 2 {x + ). (4-10) 



We look for the solution of Eq. ( |4.10| ) as a series in powers of z 2 . Note that h(x + ,x ,z) goes to 
zero proportionally to z A as z — > 0: therefore the series starts at the order z 4 and reads 



h(x + ,x , z) = z 4 h T 



(4.11) 



n=0 



Substituting Eq. ( |4.11 ) into Eq. ( [4.10 ) and solving it order-by-order in z we can express all the 
coefficients in the series in terms of the first coefficient ho(x + , x~) and in terms of t\(x~) and t 2 (x + ) 
obtaining 



h(x + ,x ,z) = 4z 2 2 ^ ^ + - h (x + ,x ) — 32 z 2 



-- z 2 d+ <9_ 
4 



24 



z 4 (d+ d. 



(d + d. 



I 2 (zy/2d + d-) 

h(x-)t 2 (x + ). 



(4.12) 



(Inverse derivatives in Eq. ( 4.12 ) are canceled by the positive powers of derivatives in the numerators 
of the appropriate terms.) 

As shown in Appendix ||, plugging Eq. (|4.12|) into Eq. (|4.9d|) one can find g(x + , x~ , z), and, using 
Eqs. Q4.91Q and (|4.9gi) , one can find f(x + , x~ , z) and f(x + ,x~, z). One can argue (see Appendix |B|) 
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that for the solution of Einstein equations to satisfy initial (h — g — f — f — before the collision) 
and boundary (h(x + ,x~,z = 0) = 0) conditions the following relation needs to be satisfied: 



(d + d^) 2 h Q (x + ,x-) = 8t 1 (x-)t 2 (x + ). 



(4.13) 



As can be seen from Eq. (|4.12|) the infinite series ( [4.11|) for h(x + , x~ , z) will then terminate at the 
order z & . As shown in Appendix |B|, the solutions for f(x + , x~ , z), f(x + , x~ , z) and g(x + , x~ , z) will 
also reduce to finite-order polynomials in z 2 . 

The only other non- vanishing coefficient in the series for h(x + ,x~ , z) in Eq. ( [4.11 ) is hi(x + ,x~) 
which is related to ho(x + ,x~) via 



hi(x + ,x ) = - d + cL. h (x + , x ). 
6 



Using Eq. ( f4.13| ) we obtain 



d + d- h\(x + , x 



(4.14) 



(4.15) 



Eqs. ( 4.13Q and (^4.15|) allow us to determine the functions ho and hi. Imposing causality by 
requiring that at time — oo, i.e. long before the collision, the shock waves are unmodified we write 



h (x + ,x ) = 8 / dx' 



X X 

J' 



dx' 1 



dx' + / dx" + t 1 (x"-)t 2 (x" + ) 



(4.16) 



and 



h 1 (x + ,x ) 



dx' / dx' + t\(x' )t2(x' + ). 



(4.17) 



In terms of h Q and hi from Eqs. ( |4.16 ) and ( ^4.17| ) we write our solution for h(x + , x , z) as 

h(x + ,x~,z) = h (x + ,x')z A + h 1 (x + ,x~)z 6 . (4.18) 

Plugging the solution ( [4 . 1 8| ) into Eq. ( |4.9d| ) we solve for g(x + , x~ , z) to obtain 

2 

g(x + ,x~,z) = -2h (x + ,x~) z 4 -2h 1 (x + ,x~) z 6 + -t 1 (x-)t 2 (x + ) z s . (4.19) 

3 

Substituting the solutions for h and g from Eqs. ( |4. 18|) and ( |4.19|) into Eq. ([4.9f|) we solve for 
f(x + , x~, z) to find 

f{x + ,x~,z) = -X^x+^-jz 4 --d 2 _h {x + ,x-)z 6 - — d^x+iX^z 8 (4.20) 

6 16 
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with Xi(x + ,x ) given by 



\i(x + ,x ) 



dx' + d_hn(x' + , x 



(4.21) 



Similarly Eq. ( |4.9g| ) yields 



f{x + ,x ,z) 



-\ 2 (x + ,x ) z 4 d 2 ,ho(x + ,x ) z t 

6 



— dlh(x + ,x )z & 
lb 



with 



(4.22) 



\ 2 {x + ,x ) 



dx' d+hn[x + ,x 



(4.23) 



Eqs. flEIgp , Q4TT9|) , ([ODD and Q02]) provide us with the solution of Eq. (Q with the causal 
initial condition requiring all these functions to go to zero at infinitely early times. 

Using Eq. fl2.3|) one can obtain the contribution of gffi to the expectation value of the energy- 
momentum tensor at the boundary of the AdS space from Eqs. ( 4.18[ ), ( |4.19| ), ( |4.2U| ) and ( [4.22D 
and Eq. ( jj : 



(T-) 



Ai(x + ,x ), (T++) 



N 2 

— %A 2 (x ,x + ), 

2 7T Z 



2 IT 2 

N 2 _ , A^ 2 
( T +-) = tt^^o^ ,x + ), (T^-) = — ^-S i:j h (x } x + ). 



2n 2 



2tt 2 



(4.24) 



Given ti(x ) and t 2 (x + ), one can use Eqs. ( |4.16|) , ( |4.21| ) and ( |4.23| ) to find h , Ai and A 2 , and then 
use Eq. ( f4.24| ) to construct the energy-momentum tensor of the gauge theory. 



5. Physical Shock Waves: Nuclear Stopping 



To understand our solution given by Eqs. flOJ), (pDD , d420|) and fl42p let us consider a specific 
example of shock waves with the boundary energy-momentum tensor given by Eq. (|3.3j ). To be able 
to better understand physical properties of the solution, let us "smear" the shock waves over some 
finite longitudinal distance. If one imagines shock waves representing a large nucleus, such nucleus 
moving in the x + -direction in a boosted ultrarelativistic frame would have a longitudinal extent 



A A 1 / 3 

a oc R— oc — — 
p+ p+ 



(5.1) 



with R the nuclear radius, A the typical transverse momentum scale in the problem (R oc A 1//3 /A 
with A the atomic number), and p + the large longitudinal momentum of the nucleus. (Here A/p + 
is the boost factor.) 
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Assuming that the nucleus has equal thickness a at all impact parameters, we replace the 
delta-function in Eq. ( |3.3|) with two theta-functions to write 



for the first nucleus and 



a 



t 2 (x+) = 2ix 2 ^9(x + )9(a-x + ) 



(5.2) 



(5.3) 



for the second one. For simplicity we assumed that the nuclei are identical and are scattering with 
equal momenta pf = p%, such that \i — //i = fi2 and a = ai = a 2 . 

Plugging Eqs. ( |5.2|) and ( |5.3|) into Eq. ([4.16 ) we immediately obtain 



h (x + ,x') = 8%(2vr 2 



x 



9{x )9{a — x )— — h 9{x — a) a [x 



x 



(x + ) 6 (a - x + ) -y- + 9(x + - a) a [x' 



(5.4) 



Eq. Q4.21D then gives 

A!(x + ,x-) = 8^(2tt 2 ) 5 



(x ) 0(a — x ) x + 0(x — a) a 



(x+) 9(a -x + ) + 9{x + - a) a 

6 



a 2 x +2 



6 + 2 



a x 



(5.5) 



while Eq. ( |4.23| ) due to the fact that nuclei are identical leads to 

A 2 (x + ,x~) = Ai(x~,x + ). (5.6) 

Eqs. (|5.4j ), (|5.5|) and (~ |5.6|) , along with Eq. (|4.24|) , give us the order /z 2 energy-momentum tensor. 
Let us study its properties. First of all, away from the light cone for x + , x~ ^> a (or in the limit of 
infinitely thin nuclei, which can be recovered by taking a — > 0) one has 



ho(x + ,x ) 



x+ ,x~ 3>a 



(27T 2 ) 2 /J 2 X + X , Ai(x + ,x ) 

A 2 (x + ,x~) 



(O 2\2 ,,2 
(2 7T ) /X 



/o 2\2 ,,2 
(2 7T J /X 



X 



(5.7) 



Substituting Eq. (|5.7| ) into Eq. (|4.24 ) one gets for the forward light-cone far away from the nuclei 



= -8 7r 2 AT 2 /i 2 x +2 , (T++) = -8 vr 2 N 2 /i 2 x' 2 , 
(T + _) = Svr^Vr 2 , (T^> = 8 vr 2 5 tj N 2 /z 2 r 2 . 



(5.8) 
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This is the same result as obtained in [67]. The energy- momentum tensor in Eq. (O) is rapidity- 
independent, as its components contain no rapidity dependence apart from the trivial factors needed 
for Lorentz-properties of the tensor. 

For two colliding nuclei of infinite transverse extent (the Bjorken case [70]) the most general 
parameterization of the rapidity-independent energy- momentum tensor is [43] 



T+- 
T = 



[<r)+p 3 (r)) 
Hr)+p 3 (r)] 



T 

- \ 2 

X 

r 



Hr)- P3 (r)]-, 



(5.9) 



where e(r), p{r) and Pzij) are the energy density, transverse pressure and longitudinal pressure 
components of the energy-momentum tensor at mid-rapidity (x 3 = 0). At x 3 = the tensor ( |5.9| ) 
looks like 



rpflV 



/e(r) \ 
p(r) 
p(r) 
y p 3 (r)J 



(5.10) 



One can easily show that conservation of energy and momentum condition 





applied to the tensor (|5.9|) gives 



de 
dr 



e + P3 



T 



(5.11) 



(5.12) 



The condition d^T^ = follows from Einstein equations if one uses gauge-gravity duality to obtain 
the energy momentum tensor. For conformal field theories the energy-momentum tensor is traceless 



rrifi 



0. 



which implies 



2p + p 3 . 



(5.13) 



(5.14) 



Eqs. ( |5.12|) and (|5.14|) relate e(r), p(r) and ^3(1") to each other, such that knowing one of these 
functions is sufficient to reconstruct the others. 
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Comparing Eq. ( |5.8| ) with Eq. ( |5.9| ) we read off the energy density at mid-rapidity 



e(r) = 4vr 2 iV c V 2 r 2 . (5.15) 

Such energy density at mid-rapidity is problematic. It grows with the proper time r. One can show 
that the requirement that the energy density of produced matter is positive-definite in any frame 
in particular demands that [54] 

e'(r) < 0. (5.16) 

Eq. ( |5.15|) obviously violates the condition (|5.16|) : this means our solution gives negative energy 
density in some frames. Such result is clearly unphysical. 

It is important to understand the origin of this negativity of the energy density. First we note 
that, as one can easily see, the energy density becomes negative in the frames with the time direction 
being close to the light cones (the shock waves). To investigate the region around the shock waves 
further, let us concentrate on the shock wave corresponding to the nucleus 1 after the collision. Let 
us study what happens to, say, the middle of the nucleus, which is located at x~ = a/2, after the 
collision. The important component of the energy-momentum tensor is (T__), since it contains 
the (large) momentum component of the nucleus along its light cone. Using Eq. ( |4.24j ) along with 
Eq. ( |5.5| ) at x~ = a/2 yields for x + ^> a (after the collision) 



(T__(x+ > a,x~ = a/2)) = N 2 C - - 4 tt 2 iV 2 /i 2 x +2 , (5.17) 



where the first term on the right is due to the original shock wave obtained by using Eqs. (|5.2|) and 
(|3) at x- = a/2. 

Eq. ( |5 . 1 7|) shows that (T__) of a nucleus becomes negative at light-cone times 



4- 1 , , 

x + — . (5.18) 

■Jji a 

Indeed zero (T_ _) would mean a complete stopping of the shock wave and the corresponding nucleus. 
We therefore conclude that negativity of energy density ( |5.15|) in fact is a signal of complete stopping 
of the colliding nuclei after the collision! 

Indeed at times x + ~ 1/yJJTa higher order corrections to the metric due to higher graviton 
exchanges would become important preventing (T__) from going negative. Nevertheless, Eq. ( [5.17 ) 



demonstrates that at rather short times x + ~ 1/ \fjla the nucleus looses the amount of energy 
comparable to its initial incoming energy, and thus is likely to stop. 

One should also point out that Eq. Q5.171 ) gives of the center of the nucleus (x~ = a/2): 

other slices of the nucleus located at different x~ would also stop, but at slightly different times x + . 



All stopping would happen at the same parametric time given by Eq. ( 5.1 
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To better understand the stopping time we use Eqs. ( |3.4j ) and ( |5.1|) to re-write Eq. (|5.18|) as 



x+ ~ Jim- ( 5 ' 19 > 

The stopping time appears to be energy-independent! It is given by the inverse of the typical 
transverse momentum scale A in the problem. It also decreases with the increasing size of the 
nucleus A. 

Let us pause to interpret the main result of this Section. It appears that two colliding ultrarela- 
tivistic shock waves would come to a complete stop shortly after the collision. One can understand 
this in terms of creation of a black hole: both shock waves carry large energy, which functions as 
a mass. There is a horizon radius in AdS space corresponding to that mass/energy. For nuclei 
of infinite transverse extent under consideration the shock waves always come closer to each other 
than the horizon radius corresponding to the energy they carry. A black hole is then formed and 
the shock waves stop completely within the black hole's horizon radius. The picture is similar to 
black hole production in collisions at transplanckian energies, which has been recently discussed in 
the literature [72-74]. 

One could picture a collision of two nuclei of finite transverse extent. If the impact parameter 
of such a collision is larger than the horizon radius, no black hole will be formed and the nuclei will 
not stop. However, in such case there will probably be no thermal matter produced in the boundary 
theory either. 

It is interesting to note that the stopping time (|5.19| ) is independent of energy: indeed, on one 



hand if one increases the momentum of the shock wave it is harder to stop it. On the other hand, 
increasing the energy of the shock waves tends to reduce the radius of the horizon, trying to make 
the shock waves stop faster. We interpret the result of Eq. (|5.19| ) as the cancellation of the two 
effects, leading to energy-independence of the stopping time. 

If the nuclei stop completely in the collision, the strong interactions between them are almost 
certain to thermalize the system. Indeed if the interactions were strong enough to stop the nu- 
clei, they should be strong enough to thermalize the resulting medium. The dynamics of such a 
rotationally-invariant thermal medium was originally described hydrodynamically by Landau in [75] 
and is commonly referred to as Landau hydrodynamics. Hence our conclusion is that modeling a 
collision of two nuclei by two physical colliding shock waves in AdS necessarily leads to complete 
nuclear stopping, and probably to thermalization of the system and the subsequent dynamics de- 
scribable by Landau hydrodynamics. 

The fact that nuclei do not stop instantaneously, but require certain (short) time ( [5.1 9| ) to stop 



avoids the standard counter-argument [89] against Landau hydrodynamics [75], which suggests that 
it would violate the uncertainty principle if the stopping was instantaneous. Hence the picture is 
intrinsically consistent. 

One may wonder whether our result of Eq. ( |5.17|) is specific to the shape of the shock waves 
we have considered in Eqs. (|5.2|) and (|5.3|). In fact our conclusion of complete stopping is valid in 
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general: as one can see from Eq. ( 4.16 ), any non- negative energy- momentum tensor of the shock 
wave, which is positive in a localized region of x~ axis, would give h ~ x + x~ at late times. 
Eqs. Q4.24D and (|579| ) would then give 



e(r) ~ p(r) ~ -p 3 (r) ~ r 2 , (5.20) 

just like in our case considered above. Hence the energy density of such system would never be 
non-negative in all frames, signaling the stopping of shock waves. Finally, Eq. ([4.21 ) would give 



Xi ~ x +2 , such that the correction to the energy-momentum tensor on the light cone would again 
be 

(T__) ~ -x +2 (5.21) 

indicating that at some large enough light-cone time x + the nucleus would run out of its momentum. 
This proves the shock waves stopping independent of the shape of the shock waves profiles. 

Let us close this Section by pointing out that the result in Eq. (|5.15|) can be easily obtained 



(at the parametric level) for infinitely thin nuclei by noticing that the diagram in Fig. |^C, which 
contributes to the metric giving the energy density in Eq. ( |5.15| ), is of the order of /i 2 . Hence the 
contribution to e at this order should be of the order of yU 2 . However, energy density has dimension 
of mass to the fourth power, while /x 2 has dimensions of mass to the sixth power. To make the 
dimensions right we use the only other dimensionful quantity in the boundary gauge theory in the 
forward light cone: the proper time r. This gives e ~ /U 2 t 2 , in agreement with Eq. ( |5.15| ). We noted 
above that expansion in graviton exchanges in the bulk is equivalent to expansion in the powers 
of \i for the energy-momentum tensor of the shock waves in Eq. (|3.3|) . The only way to make a 
dimensionless expansion parameter from fi in the boundary theory is to multiply it by r 3 . Now we 
see that for energy density the expansion parameter is in fact /it 3 , which has been noticed in [67] 
before. However, now we understand that each power of this expansion parameter corresponds to 
a graviton exchange between the boundary and the bulk. 



6. Unphysical Shock Waves: Energy Density of the Produced Medium 

In the above Section we came to the conclusion that at very strong coupling colliding nuclei com- 
pletely stop in a collision (for central collisions), forming a medium described by Landau hydro- 
dynamics. However, due to asymptotic freedom of QCD, we know that small-coupling effects play 
an important role in heavy ion collisions and in high energy collisions in general. The success of 
Color Glass Condensate [39-41] based models in describing RHIC data (see e.g. [90] and references 
therein) suggests that weakly coupled effects are present in the actual heavy ion collisions at RHIC, 
at least at very early times during and after the collision. While a comprehensive description of 
both weakly-coupled initial dynamics and strongly-coupled dynamics of the produced medium is 
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not feasible at this point, here we will suggest a model capturing some of the feature of the weakly 
coupled collisions. 

We begin by noting that, in the weak coupling limit, the colliding nuclei do not stop, as we 
observed in the previous Section for the strong coupling case. Instead the valence quarks and other 
large Bjorken-x (hard) partons are usually assumed to go through each other without deflection, 
shedding off the softer (small-x) virtual partons, which are left behind and quickly go on mass 
shell, i.e., become real [26,27,29-31,42,44,77-79]. The medium made out of these small- a; partons 
after the collision has a non-negative energy density in any frame [30,32,42,68,79]. We will then 
proceed by requiring that the energy density of the produced strongly coupled medium should also 
be non-negative. 

We want to model the heavy ion collisions by colliding two shock waves. The conclusion of the 
previous section was that any localized non-negative (T\ ) of a shock wave, such that 

oo 

J dar<Ti__(aT)) > (6.1) 

— oo 

(with an analogous condition imposed on the (T2 + + (x + )) component of the energy momentum 
tensor of the other shock wave) leads to the energy density scaling of Eq. ( |5.20| ). This violates the 
condition ( 5.16p derived in [54] and results in the negative energy density of the produced matter 



in some frames. The only way around such an unphysical conclusion appears to be to require that 

oo oo 

J dx-{T X -„{ X -)) =0, J dx + {T 2 + + {x + )) = 0. (6.2) 



Indeed the conditions ( |6.2| ) can only be satisfied in a physical world if there is no shock waves, in 
which case their energy would be zero. Such a trivial scenario is not what we have in mind. 

Instead, we propose using unphysical not positive-definite quantities for (Ti_ -(x~)) and 
(T 2 + + (x + )), which integrate out to zero satisfying Eq. (|6.2|) . Indeed such objects would be com- 
pletely unphysical, as they would contain regions of negative energy density. They can not be 
obtained from an underlying string theory either. However, we intend to use them in gauge-gravity 
duality only. On both sides of the gauge-gravity duality our non-positive energy momentum tensors 
should be regarded as external sources to the theory. The conclusion we have reached is that to 
have non-negative energy density in the forward light cone one needs unphysical negative energy 
shock waves on the light cone itself. 

Indeed our proposal of zero-energy shock waves may not be a unique way of modeling weak 
coupling effects in heavy ion collisions in the AdS / CFT framework. One may also try to construct 
a metric with the CGC-inspired energy-momentum tensor for the gauge theory at early proper time 
and evolve it in time using Einstein equations. However, constructing a metric which is a valid 
initial condition for Einstein equations at early times and accounts for perturbative features of the 
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collision appears to be difficult. If one insists on modeling the heavy ion collisions by two colliding 
shock waves, our zero-energy shock wave proposal is the only way to mimic the weak coupling effects 
at initial stages of the collision. 
Inspired by Eqs. 



HID and (3J5), which contain two factors of transverse momenta times some 
function of longitudinal coordinates and momenta, we suggest describing the energy-momentum 
tensors of the colliding nuclei by 



corresponding to 



Pi — (a")) 
<T 2 + + (* + )> 



ti(x 



N 2 
N 2 



2tt 5 



A 2 2 5\x + ) 



(6.3) 



(6.4) 



in the shock waves metric in Eq. ( [4.3| ). S'(x) denotes the derivative of a delta-function. Clearly the 
energy-momentum tensors in Eq. ( |6.3| ) satisfy Eq. ( |6.2j ). What we loose in this description is the 
relation between the transverse momentum scales A 2 and A 2 , describing the two nuclei in Eq. (|6.3| ) 
and the actual physical parameters describing the real nuclei, since our energy-momentum tensors 
in Eq. ( |6.3|) are not physical and we can not relate them to the energy-momentum tensors of the 
two nuclei. 

Before we perform any calculations we can already guess the answer using the dimensional 
analysis outlined in Sect. [|. This time each vertex in Fig. |3]C brings in a factor of A 2 and A 2 ,, 
such that the diagram is proportional to A 2 A 2 ,. Hence the resulting energy density of the boundary 
theory is proportional to e ~ A\ A 2 ,. Since the dimensions of e and A 2 A 2 , match, no powers of r are 
needed this time. Hence we conclude that the energy density of the matter produced by the two 
shock waves (|6.4|) at the lowest order in graviton exchanges is e ~ Af A 2 ,, i.e. a constant of time, as 



was suggested in [1]. e ~ A 2 Af, immediately satisfies the condition ( |5. 16| ) derived in [54]: hence the 
energy density is non-negative in any reference frame. Finally, now the graviton exchanges between 
the boundary and the bulk should correspond to powers of A\ r 2 or A 2 , r 2 in the gauge theory. Thus 
the early-time expansion for the energy density should contain powers of A\ t 2 and A 2 r 2 . Therefore, 
while we do not know how to relate A 2 and A 2 , to the physical observables, we still can systematically 
construct the dual geometry to the collision by expanding the metric in the powers of A 2 r 2 and 
AgT 2 , and hopefully would be able to arrive at the thermalization/isotropization transition in the 
Bjorken sense [43,70]. 

The actual calculations are performed easily. Plugging Eq. ( |6.4| ) into Eq. ( f4.16| ) yields 



h (x + ,x~) = 8A 2 A 2 2 6(x-) 6(x + ). 



(6.5) 
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Using Eqs. fl£T7D , CT , CT , CT) , (EH), (EH), (EH), and 
correction to the metric 



we find the second order 



(2) 

9— 



L 2 

z 

L 2 

z 



A 2 A 2 



-8 S(x~) x + 9(x+) z 4 -^ S'(x~) 9(x + ) z 6 - 5" (a; - ) 5(a; + ) z 8 



( 2 ) - ±L A2 a 2 
y++ - 2 ^i ^2 



ff? = --^A 2 A 2 

y+ 2 z 2 1 2 



<5y A? A 2 , 



y 2 U «J 



; x- 0(aT) 5(x + ) z 4 -^ 9{x~) 5\x + ) z 6 --^ S(x~) 6"(x + ) z 8 
-16 0(aT) 6(x + ) z 4 -^ 5(x~) 5(x + ) z 6 + ~ <J'(ar) 5'(x + ) z 8 



; 0(aT) 0(x + ) z 4 + - <S(aT) S(x + ) z 6 



(6.6) 



where the double prime denotes the second derivative. 



One might be concerned with the fact that order-z 4 components of g^ 2 l_ and g++ contain a 
negative contributions localized to the light cone and growing with time. They may be interpreted, 
just like Eq. ( |5.17 ) above, as a signal of stopping of the shock waves. However, our shock waves start 
out carrying negative energy and momentum densities. Hence the concept of stopping is ill-defined 
for our shock waves, because they themselves are ill-defined as physical objects, and are only used 
as some sources providing us with realistic dynamics of the produced medium in the forward light 
cone. 

Using Eqs. (EH ) an d (|5.9| ) along with Eq. (|2.3|) we can read off the energy density and the 
pressure components in the forward light cone from Eq. 



,(2) 



<r) = ^4 A 2 A 2 , 

7T 

P(r) = ^4 A 2 A 2 , 

N 2 

p 3 (r) = — f-4 A 2 A 2 . 



(6.7) 



Once again, just like in CGC [32,68] and as was obtained in [1], the strongly-coupled medium starts 
out very anisotropic, with a negative longitudinal pressure. Eqs. ( |6.6| ), combined with the lower 

allow for a systematic construction of the metric as an expansion 



order metric ( f4.3|) and Eqs. 
in graviton exchanges, to be performed elsewhere [91]. 

Negativity of the longitudinal pressure p% in Eq. ( |6.7[ ) is intimately connected with energy 
conservation. One can easily see from Eq. ( |5.12| ) that if the energy density scales as e ~ 1/r then 
Ps = 0. For the energy density which falls off slower with r, using Eq. (|5.12 ) one gets negative 
P3, in agreement with Eq. (|6.7|) . For the energy density falling off with r faster than 1/r one 
would get positive p%: however, such behavior of energy density at early time would violate energy 
conservation. The net energy of the produced medium at early times is proportional to E ~ e r. 
The energy density which increases faster than 1/r at small r would then lead to infinite energy of 
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the produced medium at very early times, violating energy conservation. Hence energy density at 
early times can not scale faster than 1/r, leading to negative or zero longitudinal pressure p 3 . 



7. The Dilaton 

Let us briefly touch upon one related topic which may become important at higher order in graviton 
exchanges. Immediately after the collision of two heavy ions the medium is not equilibrated yet. 
The magnitudes squared for the chromo-electric and chromo-magnetic fields are not equal to each 
other. This means that the expectation value of the gluonic field strength squared should not 
be zero. In fact, at weak coupling a CGC calculation at the lowest non-trivial order (order a^) 
performed along the lines of [27,43] yields 

(tri£) = -AalC F ^ I ^ [j 2 (k T r) - J 2 (k T r)] . (7.1) 

Here A and S± are the atomic number and the transverse area of the two identical colliding nuclei 
and Cp = (N 2 — 1)/2N C . To obtain Eq. ([7.1D one should substitute the gluon field from Eq. (12) 
in [27] into the Abelian part of trF 2 u and average the resulting expression in the wave functions of 
both nuclei (see also [43,92]). At early times Eq. ([7.1Q gives 



« _ 47ra ^^ ln __ ( - 2) 



where Q 2 = 4 n a 2 A/S± is the saturation scale of the nuclei which regulates the infrared divergence 
in Eq. ([7. 1|) when higher order rescatterings are included. Similar to [32,68] one may conclude that 
the scaling (|7.2|) is true for the all-order classical gluon field [26,27,29-31,42,77-79]. We get 



i 



oc-^ln-^. (7.3) 



As (tiF 2 u ) = 2(B 2 — E 2 ) with B and E the chromo-magnetic and chromo-electric fields, we conclude 
that at weak coupling the medium at the early stages of the collisions is dominated by chromo- 
electric fields (see also [93]). 

To introduce non-zero expectation value for tvF 2 in AdS one needs to include the dilaton field 
(p, as [71,82,83] 

1 < te *i> = ^ta*?tf!i>. (7.4) 



The dilaton couples to the metric through modified Einstein equations Q2.13 ) and through the 
Klein-Gordon equation ( p. 14 ). 
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The metrics of the incoming shock waves like ( |3.1|) are solutions of Einstein equations ( p. 11 ) 



with zero dilaton fields. This implies that the dilaton field is zero for a single nucleus and, using 
Eq. ( |7.4| ), (trF^ u ) = for a single nucleus as well. This agrees with the fact that at weak coupling 
(trF^ u ) = too, as the electric and magnetic fields of a single ultrarelativistic charge are equal to 
each other (the equivalent photon/gluon approximation) [24,84]. 

At strong coupling the dilaton field may become non-zero after the collision. As one can see 
from [71], an investigation of late-time dynamics of the strongly-coupled medium requires a dilaton 
field leading to non-zero (trF^ u ). Since the dilaton field of each of the shock waves is zero, we may 
only expect that the dilaton field produced in the collision would depend on energy-momentum 
tensors of both shock waves. Therefore, using the expansion in Af and A| of Sect. || one may 
expect that at the lowest non-trivial order 

(tri£> ~ A?A^ (7.5) 

corresponding to the dilaton field (see Eq. (|7.4p ) 

f(x + ,x~,z) ~ k\k\z A . (7.6) 

However, the sign of the dilaton field and (tvF^) can not be determined from these dimensional 
considerations. (In fact Eqs. ( p. 13 ) an d (|2.14j) are invariant under if — > —if and do not fix the 
sign of the dilaton field.) Instead of the logarithmic divergence Q7-3 ) of the perturbation theory, we 
anticipate the expectation value of (trF^ u ) to go to a constant at early times. 

The dilaton field from Eq. ( |7.6| ) would affect the metric only at the order Af Aj, as can be seen 
from Eq. ( [2.13 ), and therefore we could safely neglect it in the above discussion and in [1]. It would 
only enter Eq. (|2.13| ) at the same order as the four-graviton exchanges. 

We stress that while we do not know whether non-zero dilaton field would arise at higher orders 
in our expansion in graviton exchanges, if it does come in we expect it to be of the form shown in 



Eq. ( |7.6|) (at the lowest order) leading to a constant (txF^J) at early times shown in Eq. ( |7.5|) . 



8. Conclusions 

Let us summarize the main points of our paper. In Sect. [5] we have demonstrated that a collision 
of real and physical shock waves in AdS would lead to a full stopping of the shock waves. It is 
likely that a black hole would be created in AdS space immediately afterwards. For the gauge 
theory this implies that two nuclei colliding at very strong coupling would stop almost immediately 
after the collision, after a time interval of the order of t s t op ~ 1/A with A some typical transverse 
momentum of the problem. Thus for RHIC and LHC one might expect t s t op ~ 1 fm. After 
the stop, creation of the black hole in the bulk likely translates into thermalization and Landau 
hydrodynamics description of the dynamics for the created medium. This is a possible scenario 
advocated in [76]. 
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However, we consider it to be more likely that the physics of the initial stages of heavy ion colli- 
sions is weakly coupled. This point of view is supported by the many successes of CGC approaches 
to heavy ion collisions [33,34]. In the weak coupling scenario the hard (large-x) parts of the nuclear 
wave functions simply go through each other in the collisions without deflection or recoil. At the 
same time the produced thermalized medium should still be strongly-coupled. While no comprehen- 
sive single description of both the weakly-coupled early stages and the strongly-coupled medium in 
the final state exists, in Sect. ^ we constructed a model which appears to capture the main features 
of the weakly-coupled initial state by allowing the energy density to be negative (only) on the light 
cone. The corresponding shock waves carry both positive and negative energy density and are thus 
unphysical. They need to be thought about as some external sources for the gravitational field used 
in gauge-gravity duality. The energy density of the produced medium in the forward light-cone is 
non-negative: in fact we recovered our earlier result of [1] that the energy density starts out as a 
constant at early proper times. 

We have thus arrived at the following conclusion. If the coupling constant in heavy ion collisions 
is large throughout the collision this would lead to nuclear stopping followed by Landau hydrody- 
namics. In the more realistic scenario, according to present phenomenological evidence, in which 
both the nuclear wave functions and the primary particle production are weakly-coupled, Bjorken 
hydrodynamics could be still achieved if the coupling constant quickly becomes large. Indeed, as 
shown in [54], a purely strong coupling approach to the study of late-time dynamics leads to Bjorken 
hydrodynamics. However, as we argued in the paper, Bjorken hydrodynamics can not result from 
strong coupling dynamics only. 

When this paper was in the final stages of preparation, a preprint [94] was posted on the arXiv, 
where a similar conclusion about stopping of shock waves has been reached. Also, very recently a 
new version of [67] appeared on the arXiv, where the possibility of the shock wave energy-momentum 
tensor as in our Eq. ( |6.3p was briefly mentioned. 

Acknowledgments 

We would like to thank Samir Mathur for many highly educational and informative discussions. We 
thank Ulrich Heinz, Al Mueller and Krishna Rajagopal for encouraging discussions. 

This research is sponsored in part by the U.S. Department of Energy under Grant No. DE- 
FG02-05ER41377. 

A. A model for the energy- momentum tensor of an ultra-relativistic nu- 
cleus 

We begin by considering the classical electromagnetic potential of a point-like particle of charge g 
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moving at speed v along the positive z direction. In the covariant gauge it is given by 

A ± = f 0±vW2 4 = (i = 1,2). (Al) 

4 71 [\ ((1 + v )x--(l- v) x+f + (1 - v 2 ) x\] ' 

The energy-momentum tensor associated to this field is 

T, u = -F w FJ> + 1 r)^ F pa F pa , (A2) 

where r\ is the Minkowski metric in four dimensions. In the limit v — > 1, its only non- vanishing 
component is 

{ dl A.f = ^ (IZ^M 3 , (A3) 



2tt [2x- 2 + (1 



where a = g 2 /4ir. In the strict limit v = 1, T is singular at x =0. To clarify the nature of this 

singularity we consider the following integral: 



r dx € -^f(x) =(27r*)e 4 i ^ JiZ) 

J-oo (x 2 + e 2 ) 3 M J K } 2! 



2! dz 2 (z + i\e\Y 



(A4) 



z=i |e| 



where f(x) is an arbitrary analytic function that falls off at infinity rapidly enough for the previous 
integral to be well defined. From Eq. (|A3|) and Eq. ( |A4| ) we get 

T- = --^ \ (A5) 

27T37T \x~\ X]_ 

The previous equation serves as a starting point to build up a model for the energy-momentum 
tensor of an ultra-relativistic nucleus, T™ d . We envisage the nucleus as consisting of A nucleons, 
each of them containing N 2 valence gluons. Assuming an uniform distribution of nucleons inside 
the nucleus, averaging over transverse position and summing over all nucleons, we write 

<r-> = «i4-f <P« ± T__ = ^1 -4- In (*) , (A6) 

S± J ix \x \ 3b± \p J 

where S± is the nuclear transverse area, R is the nuclear radius and p is an UV cutoff introduced to 
regulate the singular behavior at x± = 0. The extra factor of N c in Eq. ( |A6| ) comes from calculating 
the color factor for non-Abelian energy-momentum tensor at the lowest order in the coupling. For 
simplicity, we consider a cylindrical nucleus such that S± ~ tt A 2 ^ R 2 ^ 3 with Rn the nucleon's 
radius. Introducing the dimensionful scale 

a2 = -Aw ln (-) ' ^ 

3 vr 3 Rjf \p J 
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we rewrite Eq. (KH) as 



(T™ d ) = N 2 C a N c 4 vr A 2 A 1/3 ^-j. (A8) 



Defining 't Hooft coupling A by Eq. (|3.6|) we rewrite Eq. (|A8|) as 



(T_ rad ) = A^ 2 A A 2 A 1/3 (A9) 



which is now ready to be cast in the form of Eq. (|3.5| ) at strong 't Hooft coupling by replacing 
A — ► vA. The singularity in 5(x~)/\x~\ can be regularized by replacing with the light-cone 

momentum of the valence gluon p + , which would reduce Eq. (|A9|) to Eq. ( |3.3|) with fi given by 
Eq. (|3.4|) . The energy- momentum tensor in Eq. (|A9|) may serve as the source for the one-graviton 
exchange shock wave metric in Eq. (|3.1|). 



B. Solution of Equations ( RET^l) 

Here we complete the solutions of the Einstein equations for the second order correction to the 
metric, g)J . Analogously to what we did for h in Eq. ( 4.12| ), we write the unknown functions g, f 



and / in Eq. ( |4.8| ) in the form of a power series in z 2 , starting at order 



DC 

g(x + ,x~,z) = z 4 ^2g n (x + ,x~) z 2n 

n=0 

oo 

„2n 



[X'X Z 



f(x + ,X ,Z) = Z 4 J^/„( 

n=0 

oo 

f(x + ,x~,z) = z 4 ^/„(x + ,x")z 2ra . (Bl) 



n=0 

Inserting these expansions into Eqs. ( f4.9e|) , ([4.91]) and ( |4.9g| ) we find 



2 

g n + 2h n = -6 n , 2 t l (x-)t 2 (x + ), (B2) 
f n ,x- + K :X + = -— S n>2 h(x~) t' 2 (x + ) , (B3) 

fn,x+ + K,x- =-Jo Sn ' 2 ^( x ~) ^(> + ) , (B4) 

which straightforwardly lead to the following relation between the metric coefficients 

2 

g(x + , x~, z) = -2 h(x + , x~, z) + - z 8 ti(aT) t 2 (x + ) , (B5) 

3 

f(x + ,x-,z) = -^(h(x + ,x-,z) + ^zH 1 (x-)t 2 (x + yj , (B6) 



f{x + ,x ,z) 



_|t (h(x + ,x-,z) + ±z 8 t 1 (x-)t 2 (x + )\ . (B7) 
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To complete the solution of Einstein equations Eqs. ( [4.9| ) we insert the relations in Eqs. 
B7| ) into the remaining Einstein equations, Eqs. (|4.9aLfOk| ), which we have not used yet. We find 
that, while Eqs. ( |4.9a| ) and ( |4.9bj) are trivially satisfied, Eq. ( |4.9c| ) provides the following equation 
for h(x + , x~ , z): 



3h z — z h zz + 2 z h x -\ 



z t\\x )t 2 (x 



which can be easily reduced to Eq. (|4. 10| ) with the general solution given in Eq. ( |4.12j ). Finally, 
in order to complete the solution of the Einstein equations ( |4.9|) we must determine the unknown 
coefficient ho(x + ,x~) in Eq. ([4.12|) . We do so by imposing causality of the solution which, as 
discussed in Sect. |J implies that all the second order corrections to the metric in Eq. ( }4.8|) must 
be zero before the collision. Given the relation between the metric coefficients in Eqs. ( |E5| - |F3T| ) , it 
is sufficient to require causality of h(x + ,x~,z). We first note that by requiring that 



(d+d-) 2 h {x + ,x-) = 8t 1 {x-)t 2 {x + ) 



(B9) 



we satisfy the initial (h = before the collision) and boundary (h(x + ,x~,z = 0) = 0) conditions 
on the solution of Eq. (|B8|). Eq. ( P9| ) is exactly Eq. ( f4.13|) above. This condition makes the infinite 
series in Eqs. (|4 .12 ) and ( |B"T|) terminate, as previously announced, and satisfies the causal initial 
conditions. As the solution of Eq. ( |B"8| ) for a given set of initial and boundary conditions is unique, 
we conclude that the condition ( |BTS| ) singles out the only causal solution of Einstein equations. 
(This implies that solutions obtained from Eq. ( 4.12| ) without terminating the infinite series are 
not causal.) This proves that Eq. ([B9|) gives us the only solution of Eq. (|B8| ) satisfying the needed 
initial and boundary conditions. The complete solution of Einstein equations ( }4.9|) is given by Eqs. 
( HI ), Q , (j]P and jgp . 
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